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Abstract. The structure of the near-polar accretion spots on accreting 
magnetic white dwarfs has been studied theoretically and observationally 
in numerous papers over the last decade. Detailed treatments are avail- 
able for the regime of low mass flux, usually termed the bombardment 
case, and for higher mass fluxes which create a strong shock standing 
above the photosphere of the white dwarf. No general treatment is so far 
available for the case of shocks buried deep in the photosphere. I review 
the theoretical foundations, present some applications of theory, and dis- 
cuss in short the open questions which still need to be addressed. 



1. Introduction 

Observational evidence indicates that the near-polar accretion regions on the 
white dwarfs at least in some magnetic CVs (mCVs) are highly structured, 
displaying rapid fluctuations of the mass flux both in the spatial coordiates and 
in time. Since information of their structure is lacking, a look at well-studied 
cases of magnetic accretion in the solar system may be inspiring: (i) Jupiter's 
auroral oval displays substantial structure in its UV emission as observed with 
the Hubble Space Telescope, and (ii) the Earth's aurora shows narrow striations 
which extend along field lines and indicate different densities in adjacent flux 
tubes. Although the physics is different, similar structures may exist in mCVs. 

Present models of mCV accretion spots are not sufficiently sophisticated to 
account for such structure. They assume a spatially constant mass flux across 
the column and are respresentative rather of an individual subcolumn than of 
the accretion region as a whole. The radiative interaction between individual 
subcolumns represents a complication which is usually avoided. This neglect 
may be justified in the optically thin case, but represents a severe restriction if 
the column is optically thick. In this review, I will summarize the efforts made so 
far in understanding the properties of accretion columns with special emphasis 
on the solution of the coupled hydrodynamic and radiative transfer equations. 



2. The Accretion Scenario 

Figure la schematically shows the pole cap of the magnetic accreting white 
dwarf in a polar. The accretion column is fed by a mass flux m = poVo (cgs unit 
gcm~^s~^) with po the density and Vq the value of the frre-fall velocity vg at the 
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Figure 1. (a) Left: Schematic diagram of a (tall) accretion region 
on the white dwarf (from Gansicke, priv. comm.). (b) Right: Overall 
spectral energy distribution of the prototype polar AM Herculis at 
orbital phase (p ~ 0.4. The dots represent the UBVRIJHK fluxes, 
the soft X-ray spectrum is from the April 1991 ROSAT observation, 
and the lUE UV and the HEAO-1 hard X-ray spectrum are shown 
by solid lines. The dashed line represents the unabsorbed soft X-ray 
(quasi-)blackbody spectrum. 

shock front. The size of the accretion spot depends on the coupling processes 
in the Li-point and in the magnetosphere of the white dwarf, while the stand- 
off distance of the shock front, /i ^ | i^o^cooi with tcooi the cooling time scale, 
depends on the emission process, i.e. on the electron temperature Tg, and the 
magnetic field strength B. Depending on the parameters, the emission region 
in mCVs is either tall (as shown in Fig. 1) or pill-box shaped (i.e. smaller than 
wide) and the dominant emission process is either hard X-ray bremsstrahlung 
or IR/optical cyclotron emission. Reprocessing in the atmosphere of the white 
dwarf produces, in addition, soft X-ray emission from the immediate surrounding 
of the impact region and UV emission from a wider irradiated area. 

The fact that some polars are dominated by soft X-ray emission has been 
discussed in the literature for some 20 years and dubbed the "soft X-ray prob- 
lem". In the original model by Lamb & Masters (1979) and King & Lasota 
(1979), the white dwarf was assumed to intercept ^ 50% of the emission, i.e. 
the non-reflected fraction of the downward flux, and re-radiate it as soft X-rays. 
Kuipers & Pringle (1982) effectively solved the mystery of the intense soft X-ray 
emission by noting that dense blobs of matter would penetrate the photosphere 
and the subsequent optically thick radiative transfer would reprocess all emis- 
sions into quasi-blackbody soft X-rays (see alo Frank et al. 1988). 

As an example. Figure lb displays the overall spectral energy distribution 
of the prototype polar AM Herculis at an orbital phase at which near-maximum 
cyclotron and X-ray emission is seen. The spectral energy distribution nicely 
displays the expected components: (i) cyclotron emission at infrared and opti- 
cal wavelengths; (ii) the heated polar cap of the irradiated white dwarf in the 
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UV longward of the Lyman edge; (iii) an intense soft X-ray component; and 
(iv) a hard X-ray bremsstrahlung (-|- emission hne) component. The companion 
star and the unheated fraction of the white dwarf are comparatively faint Hght 
sources. All accretion-induced emissions fluctuate in time. E.g., the soft X-ray 
flux of AM Herculis in its high state displays statistically relevant variations 
down to a time scale of ~ 0.2 s, as shown by ROSAT PSPC observations. This 
time scale corresponds to density fluctuations in the free-falling pre-shock ac- 
cretion flow on length scales as short as 10^ cm. An earlier interpretation of the 
EXOSAT soft X-ray light curve of AM Her (Hamcury & King 1988) suggested 
already that many such intcpcndcnt subcolumns are present at any time. For 
simplicity, most models disregard the inhomogeneities in the pre-shock flow in 
space and time. 

For sufficiently high mass flux, a standing shock separates the infalling 
matter from the shocked plasma which slowly cools and finally settles onto the 
white dwarf. An analytic solution for the structure of an optically thin postshock 
flow cooling by bremsstrahlung was given by Aizu (1973). If the flow is optically 
thick, which is always the case for cyclotron emission in the fundamental, the 
solution becomes much more complicated because the radiative transfer has to 
be treated simultaneously with the hydrodynamics (Woelk &; Beuermann 1996). 



3. Bombardment solution 

In a hydrodynamic shock, the ions are stopped by Coulomb collisions and the 
shock thickness equals about one ion mean free path. A special situation arises if 
cooling is so strong that the randomized energy can be radiated by the electrons 
within that scale length, i.e., effectively "on the spot". In this case, an extended 
post-shock cooling flow does not exist and the concept of a discontinuity sepa- 
rating two regimes of the flow is no longer appropriate. The whole flow has now 
collapsed to a width which corresponds to the shock thickness and effectively 
can be described as a hot corona forming part of the white dwarf atmosphere. 
This mode is referred to as the bombardment solution. 

The required efficient cooling for this case can be provided by cyclotron ra- 
diation in a strong magnetic field which causes the electron temperature to stay 
far below the Rankine-Hugoniot temperature Tgh of the shock solution. Orig- 
inally suggested by Kuijpers & Pringle (1982), this case was rigorously solved 
by Woelk & Beuermann (1992). They treated the ion collisional energy loss by 
means of statistical plasma physics and simultaneously accounted for the losses 
of the electrons by a fully frequency-dependent and angle-dependent radiative 
transfer in a plane-parallel heated atmosphere. The solution predicts a peak elec- 
tron temperature which increases as kT^ oc mP'^"^ (Woelk &i Beuermann 1993, 
see also Fischer Sz Beuermann 2001) and equals IkeV for rh ~ 0.01 gcm~^s~^, 
B ~ 30 MG, and a white dwarf mass of 0.6 Mq. Such low temperatures arc ob- 
served in AM Her stars with very low accretion rates (Schwope this conference) 
and suggest that the predicted temperatures are roughly correct. The solution 
becomes invalid when Te approaches the Rankine-Hugoniot temperature with 
increasing m. 
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4. Radiation-hydrodynamics 

For rh ^ 0.1gcm~^s~^, the post-shock plasma flow and the radiative transfer 
have to be treated simultaneously. The only feasible approach is a fluid de- 
scription which does not allow to resolve the shock front, however. There is no 
simple connection between such an approach and the kinetic description of the 
bombardment solution. Nevertheless, we find that the fluid approach for low m 
recovers some important apects of the bombardment case. 

4.1. The Equations 

The general case assumes a two-fluid plasma and allows for different temper- 
atures Ti^e and bulk velocities fi^g of ions and electrons. At the high densities 
in mCV shocks. Coulomb collisions will assure Vi = = v, except perhaps 
within the (unresolved) shock front. A two-fluid treatment may then proceed 
with a single equation of motion which contains the total pressure, i.e. the sum 
of ion pressure Pi and electron pressure Pg- The skalar pressure and the tem- 
perature are related by Pi^c = n-[^ckTi.c (with n-i and n-e the particle densities of 
ions and electrons), which assumes that the particle distributions are isotropic 
and Maxwellian. Cyclotron losses could, in fact, produce an anisotropic electron 
distribution function, but the optical depth in the fundamental harmonic is so 
large that radiative transfer causes the temperatures perpendicular and parallel 
to the magnetic field to agree closely, except very near the surface of the emis- 
sion volume (Zhelezniakov 1983). For simplicity, we assume a hydrogen plasma 
in which n[ = and the mean molecular weight is = 1/2. 

The hydrodynamic equations include the continuity equation, the equation 
of motion, and the energy equations for ions and electrons, which describe the 
conservation of mass, momentum, and energy, respectively: 

f+p(V-^l = (1) 

5 + -V(Pi + Pe) = 5 + Aad (2) 

at p 

dPi P\dp , ^. , ^ , . 

^ + 7-^ = (7 - 1) [a., + Ai - V . (,-+ F„i)l . (4) 

at p at L J 

Here, ji = -§i+v-V,p = m^rii is the mass density with mu the unit mass 
(~ proton mass), and 7 = 5/3 is the adiabatic index. On the right hand sides, 
the equation of motion contains the gravity g and the volume force exerted by 
absorbed and scattered radiation, /rad- The latter becomes important when the 
accretion rate approaches the Eddington limit. In the energy equations, Aci 
describes the energy exchange between ions and electrons by Coulomb collisons 
and Ah refers to electron heating by any other process, with the energy provided 
again by the ions. Additional energy gains of the electrons are described by 
the divergence of the heat conduction fiow q and of the radiative fiux Prad- 
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Calculation of the radiative flux (Eq. 5) requires solution of the radiative transfer 
equation (6) for the intensity Ii,{n,u) as a function of the direction n and the 
frequency where a and % are the scattering and absorption coefficients, Jy is 
the angle-averaged intensity and the use of Bi, for the source function indicates 
the assumption of LTE: 



Unsurprisingly, Eqs. (l)-(6) have not been solved in all their beauty, but only 
after the introduction of some rather severe simplifications. 

4.2. Simplified Geometry 

The correct funnel geometry (see Canalle, this conference) is most complicated 
so solve, the spherical or the linear geometry represent the simpler cases of 
which the latter is an acceptable approximation only for shock heights /ish ^ 
i?wd- Realistic columns are always limited in the direction perpendicular to 
the flow. This limitation is of no concern when the radiative losses can be 
treated as optically thin emission, but in the optically thick case it requires a 
multi-dimensional radiative transfer (i.e., with the direction of VTe position- 
dependent). Comparatively simple cases are (i) a shallow infinite layer with 
the radiative flux vector constantly pointing in the radial direction (Woelk & 
Beuermann 1996, Fischer &: Beuermann 2001) or (ii) a tall column with the 
radiative flux vector pointing only sideways (Cropper et al. 1999). Naturally, 
gravity effects increase with column height. 

4.3. Electron Heating by Non-Coulomb Processes 

In Eqs. (3) and (4), the quantity Ah represents any heating process other than 
Coulomb collisions. Non-Coulomb heating of electrons is known to occur in in- 
terplanetary shock fronts, in particular, the bow shock of the Earth (e.g. Stone 
& Tsurutani 1985). In these collisionless shocks (with a typical particle den- 
sity of n ~ 10cm~^) characteristic processes are the reflection of ions, charge 
separation between electrons and ions, and the incidence of various types of 
plasma instabilities which lead to rapid electron heating. Plasma instabilities 
may also play a role in heating the electrons in supernova shocks (Lesch 1990). 
Whether such processes occur in mCV shocks with v \\ B, high Mach numbers, 
and high particle densities (exceeding those in interplanetary space by factors 
of 10^^ to 10^^) needs to be studied. Some authors (e.g., Saxton & Wu 2001) 
assume the presence of plasma heating in the shock front and consider Tc im- 
mediately behind the shock a free parameter, possibly ranging up to Ti. Others, 
like Woelk Sz Beuermann 1992, 1996) and Fischer & Beuermann (2001) have 
disregarded plasma heating and considered Coulomb collisions only (A^ = 0). 
In their solutions, the ion temperature displays a (quasi-) discontinuity at the 
shock, while the electron temperature increases only on a length scale which 
equals the Coulomb time scale times the flow speed. Clearly, the question of 
plasma heating in the shock front needs further consideration. 
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4.4. One-fluid Plasma 

Electron and ion temperatures will be similar if the time scale for energy ex- 
change between ions and electrons by Coulomb collisions is short compared with 
the cooling time scale. A one-fluid plasma is a reasonable approximation for high 
m and low B, while a two-fluid approach is required for low m and high B. The 
dividing line between the two regimes is approximately located at mB^^'^ ~ 0.1, 
where rh is in gcm~^s~^ and B-j is in units of 10^ G (see Fig. 4 below). Hence, 
Te ~ Ti for mSf^-^ > 0.1 and < Ti for mSf ^'^ < 0.1. 

4.5. Cooling functions 

The free-free optical depth of the post-shock flow parallel to -u is rg ~ 0.1 . . . 1.0. 
Hence, a flow cooling by bremsstrahlung only is optically thin to a flrst approxi- 
mation with V-Frad = 1-4 lO^^'^n^Tc^^^ crgcm~^s~^. In this case, the divergence 
term is a locally defined quantity which obviates the need for solving the radia- 
tive transfer equation. Most authors have simplified the radiative transfer by 
replacing V • Fj-ad by a generalized cooling function expressed as the sum of terms 
of the form p'-^T^ oc p'^^^P^ with a, b fixed for each term. Bremsstrahlung is rep- 
resented by a = 2,6 = 0.5, cyclotron radiation by a = 0.15,6 = 2.5 (Saxton & 
Wu 2001, and references therein). The concept of a cooling function is applica- 
ble in approximately the same limit as the one-fluid plasma, mBj'^'^ ^0.1 (see 
Fig. 4 below) . This limitation applies to radiative transfer through the shock 
front. For tall columns which lose energy preferentially through their sides, 
cooling functions are applicable also for lower values of rhB^^'^ 

4.6. Energy Transport by Conduction 

5 /2 

The heat conduction flow, q oc UeTc VTe, is of importance only in regions of 
high particle density, high electron temperature, and/or a large gradient of the 
latter, i.e. (i) in the transition between the cooling flow and the stellar atmo- 
sphere and (ii) in the shock front with the adjoining pre-shock flow (Imamura et 
al. 1987). Most attempts to solve Eqs. (l)-(6) proceed with q = 0. Wu (2000) 
has summarized the effects of heat conduction in mCVs. 

4.7. Atomic Line Emission 

The post-shock region is optically thick for the stronger atomic emission lines, 
complicating their calculation. One may adapt a procedure known from stellar 
atmosphere research which involves two steps: (i) the temperature structure is 
calculated assuming continuum absorption and emission only; (ii) the emerging 
emission line fluxes are calculated using this temperature structure. Tennant et 
al. (1998) and Wu et al. (2001) have used a cooling function for step (i) and a 
simplified line transfer for step (ii). An accurate calculation of the line profiles 
and line intensities requires substantial effort. 

5. Oscillatory Solutions and Linear Stability Analysis 

Time-dependent solutions of the hydrodynamic equations can be used to study 
the stability of stationary solutions. The added complication of the time depen- 
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dence is compensated for by the simplicity of the coolong function. The stabiUty 
analysis proceeds by, e.g., (i) assuming small perturbations of the position of the 
shock front and of the plasma parameters, (ii) linearization, and (iii) analysis of 
the growth rates of the harmonics of the oscillations (Langer et al. 1982, Saxton 
& Wu 2001, Wu 2000 and references therein). The possible frequencies are the 
equivalent to the harmonic oscillations of a pipe which is open at one end. 

Saxton & Wu (2001) confirm earlier findings that shocks cooling by brcms- 
strahlung only are unstable against oscillations, except possibly in the funda- 
mental mode. Shocks dominated by cyclotron cooling, on the other hand, are 
stable in the fundamental and all reasonably occurring higher oscillation modes. 
This suggests that stationarity is a reasonable model assumption if cyclotron 
emission is the dominant cooling agent. The low-amplitude oscillations with 
periods of ~ 1 — 3 s seen in several mCVs (e.g. Larsson 1982) probably originate 
from regions within the accretion spot which arc dominated by bremsstrahlung. 

Solutions of Eqs. (l)-(6) with a time-dependent mass flux m incident on the 
white dwarf do not exist to my knowledge. 



6. Stationary Solutions with Cooling Function 

With V • -Frad = Ac oc p3.-bpb ^^j. g^j^ q£ g^^]^ terms), only the hydrodynamic 
equations have to be solved. In the stationary case for linear or spherical geom- 
etry, Eqs. (l)-(4) reduce to a system of coupled ordinary (though non-linear) 
equations, in which all plasma parameters depend only on the radial coordiate. 

6.1. One-fluid Plasma 

For a one-fluid plasma, P = Pe + Pi and Eqs. (3) and (4) may be combined to a 
single energy equation. The linear case is simplest and applies if the shock height 

is hsh < Rwd- With /rad = 0, Ah = 0, g*= 0, 7 = 5/3, g ~ const, z = r - i?wd, 
and f < 0, the hydrodynamic equations reduce to 

A„„) = o (7) 

^cfo^l(lP__^ (8) 
dz p dz 

V—--y-^- l\ (9) 

dz 3 p dz 3 ^' 

For the above cooling function Ac, a closed- integral solution has been given Wu 
et al. (1994). Because of the one-fluid assumption, the solution is limited to the 
high rh, low B part of the m-B plane. In the spherical case, the mass flux pv 
in Eq. (7) and g in Eq. (8) are r-dependent. The linear and spherical cases can 
be solved analytically if cooling is by bremsstrahlung only (Aizu 1973). 

A potentially important case is that of a narrow, tall column (width ^ 
height <C Rwd) which emits cyclotron radiation primarily sideways (_L B). The 
radiative flux vector is then approximately in the x, y-plane which effectively 
decouples hydrodynamics and radiative transfer. With certain simplifying as- 
sumptions on the transverse radiative transfer, the radiative losses can be ap- 
proximated by a cooling function (Wu et al. 1994, Cropper et al. 1999). 
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6.2. Two-fluid Plasma 

For strong cyclotron cooling {rriBj'^'^ <C 0.1), the one-fluid approximation be- 
comes invalid. Electron and ion temperatures differ, and Eq. (9) has to be re- 
placed by 

dPi 5 Pi dp 2 , . . , 

Although the temperatures deduced this way will be more reliable, the question 
of electron heating in the shock (Ah) implies a remaining uncertainty. Temper- 
ature profiles of ions and electrons with different assumptions on Ah (and po) 
are depicted, e.g., in Saxton &: Wu (2001, their Fig. 3). For strong heating, 
the profiles are characterized by a maximum at the shock and a rapid decrease 
thereafter, while weak shock heating causes the electron temperature to assume 
a maximum further downstream, approaching the case of pure Coulomb heating. 



7. Stationary Plane-parallel Two-fluid Radiation-hydrodynamics 

The only truely radiation-hydrodynamic solution is that of Woelk &; Beuer- 
mann (1996). Fischer & Beuermann (2001) present further results including a 
parametrization of the temperature profiles. 

Their approach involves the solution of the four hydrodynamic equations, 
Eqs. (7), (8), (10), and (4) with Ah = and g = 

dz 3 p dz 3 V dz J 

simultaneous with the radiative transfer equations (5) and (6). The solution of 
the latter was developed from a Feautrier code for plane-parallel stellar atmo- 
spheres, modified to include the anisotropy of cyclotron absorption and emission. 
The radiative transfer is fully frequency and angle dependent. The solution is 
energy conserving and ensures that the radiative flux (Eq. 6) exactly equals the 
accretion energy ^mv"^. However, because of the plane-parallel ansatz the ra- 
diative flux is directed radially. In his review, Wu (2000) incorrectly states that 
the diffusion approximation was used and that energy is not conserved. 

The solution relates Te in a unique way to m, B, and Mwd as the physical 
input parameters. Fig. 2 shows examples of the proflles of Te{z), Ti(z), and v{z) 
for B = 30 MG and m = 0.01,0.1 and 1.0gcm~^s~^. Electron and ion tempera- 
tures are found to differ substantially near the shock front and equilibrate only 
further downstream. All flows are optically thick in the cyclotron fundamen- 
tal which requires a negative electron temperature gradient to allow the energy 
to escape through the shock front. This is a consequence of the assumption 
of an infinite layer and may no longer hold true for columns of finite width. 
The effect of plama instabilities heating the electrons in the shock has not so 
far been considered by us and one can only speculate as to the effect they may 
have in the presence of optically thick radiative transfer. If the radiation escapes 
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Figure 2. Profiles of T[, T^, and v for B = 30 MG, and m = 
0.01, 0.1, and l.Ogcm^^s^^ (from Fischer & Beuermann 2001). 
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Figure 3. (a) Left: Peak electron temperature in units of the ion 
temperature at the shock for various values of the field strength vs. 
mBj'^'^. (h) Right: Normalized shock height hg^Bf^ vs. mBj^'^. 



through the shock front, radiative transfer is likely to counteract plasma heating 
and will tend to depress Tg. I consider it unlikely that Tg at the shock can be 
anything near the Rankine-Hugoniot temperature for high B and low m. The 
low temperatures of ~ 1 keV observed from some AM Her stars (e.g. Schwope, 
this conference) seem to support this conclusion, but it is fair to state that the 
importance of plasma instabilities in mCV shocks is far from understood and 
further study is needed. 

In Fischer & Beuermann (2001), we have parameterized the profiles Te{z) 
by their peak value Tmax and the shock height h^h- Fig. 3 shows Tmax normal- 
ized to the ion temperature at the shock which initially contains all the energy 
(left panel) and a normalized version of the shock height (right panel). Both 
representations unify the results for a large range of model parameters ( "data" 
points) into a unique dependence on the quantity rriBj'^'^ ("data" points and 
solid curves for 0.6 Mq, dotted curves for IMq). It is noteworthy that the de- 
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creasing peak electron temperature for low rriB^"^'^ merges into the results for 
the bombardment solution (straight lines). In the same limit, the normalized 
shock height /igh^r'^ (in cm) becomes constant. From Fig. 3 (right panel) one 
finds a limiting shock height for the 0.6 white dwarf and for low rriB^^'^ of 
hsh ^ l.giO^Bf^ ^ i.e. 10^ cm for 30 MG and as low as 5 10^ cm for 100 MG. 
Since bremsstrahlung-dominated shocks at the same B but higher rh will be still 
lower, the assumption of plane-parallel or pillbox-shaped shocks is likely to be 
good for the higer field strengths. 

In an attempt to generalize the results for a finite column width D, we 
have calculated the additional energy losses by optically thick radiation leaving 
their sides using a ray tracing technique (Fischer & Beuermann 2001). For 
a given m, these added losses cause the equilibrium temperature to be lower 
than in an infinite layer. A correct treatment of the problem would require 
a truly 2-D radiative transfer which presents a substantial complication and 
introduces the with D as an additional model parameter besides m. No such 
calculations are available in the literature. Fischer Sz Beuermann suggest a 
simple correction which ensures energy conservation: the electron temperature 
derived for for an infinite layer approximately applies to the higher value 
mD = rmoo for finite D, where r = jO/rhVo is the ratio of the total radiative loss 
from the finite region over the accretion energy (both per cm^). They provide 
a calibration of r as a function of /igh/D mB~^-^. Finally, note that in all 
models which approximate the optically thick radiative transfer by an emergent 
Planck spectrum cut off at some critical frequency v* energy is, by definition, 
only approximately preserved. 

8. Applications 

8.1. Free-standing vs. Buried Shocks 

The ram pressure Pram = Po^o the infalling matter suppresses the bottom of 
the cooling flow into the atmosphere of the white dwarf. If also the shock front 
is located within the atmosphere, the shock is considered submerged or buried 
(Frank et al. 1988). Fig. 4 compares hgh vs. m (from Fig. 3, right panel) to the 
depression produced by Pram- one finds that bremsstrahlung-dominated shocks 
are buried for rh ^ lOgcm"^, while cyclotron-dominated shocks in high-field 
polars may be buried for much lower m. To be sure, the chance for the upward 
directed radiation to escape will depend on the details of the geometry. 

8.2. Emerging spectra as a function of m 

We have employed a heuristic model to estimate the effect of buried shocks on 

the emerging spectra. Specifically, we use a field strength of P = 14 MG (appro- 
priate for AM Herculis) and consider radiation emerging at an angle against the 
radial direction of 60°. We consider subcolumns with rh values ranging from 10~^ 
to lOOgcm^^s^^, disregard the radiative interaction of the subsolumns, adopt a 
warm absorber which leaves soft X-rays largely unaffected, and approximately 
account for photoelectric absorption in the surrounding atmosphere using a ray 
tracing technique. Radiation intercepted by the white dwarf atmosphere is as- 
sumed to be reprocessed into UV radiation and soft X-rays with a temperature 
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Figure 4. Shock height vs. mass flux m for a white dwarf mass of 
0.6 Mq and five values of the polar field strength. Also shown is the 
depression of the bottom of the post-shock flow below the photosphere 
(photo) and the level in the undisturbed atmosphere which corresponds 
to optical depths r = 1 for 1 keV and 10 keV photons (Beuermann &: 
Fischer in preparation). 




Frequency (Hz) 

Figure 5. Representative emerging spectra for B = 14 MG, = 60° 
against the field direction, and mass fluxes m = 10~^, 10~^, 10~^, 1, 10, 
and 100gcm~^s~^, each for an emitting area of 10''^^ cm^ at 100 pc dis- 
tance. Dashed curves are the emitted spectra, solid curves the emerging 
spectra. The blackbody assumption is a crude approximation, in par- 
ticular for the heated white dwarf at low rh (Beuermann &; Fischer in 
preparation). 
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Figure 6. (a) Left: Synthesis of optical/IR cyclotron fluxes of AM 
Her in its high state (Priedhorsky et al. 1978). (b) Right: Synthesis of 
ROSAT and HEAO-1 X-ray spectrum (Gansicke et al. 1995, Rothschild 
et al. 1981) (from Beuermann &; Fischer in preparation). 



at the base of the individual subcolumns and their surroundings, "which depends 
on the energy flux incident on the "\vhitc d"warf. For simplicity, wc assume black- 
bodies for the reprocessed flux. This assumption may be acceptable for the hot 
sections of the atmosphere, but the cooler UV emitting parts of the heated pole 
cap "would correctly sho"w the Lyman edge in absorption (Gansicke et al. 1998, 
Fig. 12). In Fig. 5, the original emission and the reprocessed spectra are shown 
as dashed and solid curves, respectively. For low rh, the shock stands above the 
atmosphere and the emitted radiation escapes more or less freely, for larger m 
the originally emitted radiation is increasingly absorbed and reprocessed into 
soft X-rays, until at m = 100gcm~^s~^ the emission region is hidden behind an 
average column density of ~ 10^^ gcm~^ and is completely reprocessed. 

8.3. Spectral synthesis of AM Herculis 

The characteristic shapes of the spectra in Fig. 5 depend on rh and B. For a given 
field strength, such set of spectra can be used to synthesize the observed spectral 
energy distribution and thus yield the generating m-distribution. Fig. 6 shows 
the result for AM Herculis, based on the cyclotron and X-ray spectral energy 
distributions observed near the respective orbital maxima. The cyclotron flux 
(Priedhorsky et al. 1978) is modelled by contributions with 100gcm~^s~^ ( 
curve a) to rh = 10~^gcm~^s~^ (curve f), with very small contributions for 
rh ^ 3gcm~^s^^. The hard X-ray spectrum, on the other hand, originates 
mostly from mass fluxes around 1 gcm^^s^^, and the large soft X-ray flux is 
produced by rh ;^ 10 gcm~^s~^. Hence, the origin of the individual contributions 
is to a sufficient extent mutually exclusive. It is comforting and adds to the 
credibility of the result that the contributions for rh around 1 gcm^^s^^ derived 
from the cyclotron and X-ray fluxes are consistent with each other. The same 
analysis has been performed for the low state of AM Herculis, using the infrared 
cyclotron spectrum of Bailey et al. (1991) and the low-state ROSAT PSPC 
spectrum of September 1991. 



Radiation-hydrodynamic Models 



13 



10" 























I 





lO'' lo' 10^ lO' 10° lo' 10^ lo' 

m (g s cm ) 



Figure 7. Contributions AM to the accretion rate M from mass 
fluxes m between 10~^ and 100gcm~^s~^ for AM Herculis in its high 
state (filled circles) and its low state (open circles) (Beuermann & Fis- 
cher in preparation). 



8.4. Distribution of Mass Fluxes in the Spot on AM Herculis 

Figure 7 depicts the contributions AM to the total time-averaged M in AM Her 
in its high state as derived from the spectra in Fig. 6. Also shown is the result 
of a corresponding analysis of the low state. While the errors in the derived 
values of AM are much larger in the low state, a plausible picture emerges. 
The high state is characterized by a large contribution from high m which are 
largely missing in the low state, while a drizzle of low m seems to be present 
at all times (including most low states observed so far). Quite possibly, these 
differences exist already in the flow passing through the nozzle in Li. Note that 
"100gcm~^s~^" is actually a synonym for dense blobs since we are not able 
to distinguish between more or less dense blobs. The fate of blobs of different 
density is not considered here as long as the shocks are likely to be buried 
(compare Frank et al. 1988). 

8.5. Variation of Hard vs. Soft X-rax Fluxes 

It is long known that soft X-ray emission does not prevail in all AM Her stars, 
a notable exception being some systems with a low polar field strength (Beuer- 
mann Sz Schwope 1994, Ramsay et al. 1994, Ramsay this conference). One possi- 
ble explanation is the shift of the whole m-distribution towards smaller m, which 
is expected to occur in systems with lower field strengths (Beuermann 1998). 



9. Outlook 

An attempt to model the entire internal structure of the emission spots in po- 
lars is probably futile. Some important aspects of the relevant physics need 
additional attention, however. Among them are the possible heating of elec- 
trons by plasma instabilities in mCV shocks, the effect of the optically thick 
radiative transfer on T^, and the further development of two-fluid models. Pio- 
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neering studies of two-dimensional radiative transfer may be enlightening, too. 
On the observational side, studies high time resolution with the large ground- 
based and space-borne telescopes will undoubtedly provide a rapidly increasing 
insight into the dynamics and the spatial structure of the cyclotron as well as 
the X-ray emission regions. 
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